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ABSTRACT 


A  method  of  simulating  a  thin  plasma  sheath  of  the  type  encountered  in  some 
aerospace  applications  by  using  a  single  plane  of  equally  spaced  wires  is 
described.  It  is  shown  that  the  surface  impedance  of  such  a  periodic  struc¬ 
ture  can  be  made  equal  to  that  of  a  thin  plasma  sheath.  In  this  treatment, 
the  respective  surface  impedances  are  related  through  three  dimensionless 
parameters.  A  simple  technique  for  obtaining  high  losses  is  also  discussed, 
whereby  collision  frequencies  in  the  neighborhood  of  the  operating  frequency 
are  realizable. 

Experimental  data,  including  both  the  lossy  and  lossless  cases,  are  compared 
with  the  theoretical  model  in  the  X-band  region. 
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I.  INTRODUCTION 


The  simulation  of  the  electromagnetic  properties  of  a  plasma  sheath,  such  as 
those  encountered  with  a  slender  body  re-entry  vehicle,  would  be  useful  in 
laboratory  antenna  pattern  measurements.  The  main  objective  of  this  report 
is  to  develop  a  technique  of  simulating  the  reflection  and  transmission  charac¬ 
teristics  of  a  very  thin  sheath.  The  sheath  will  be  assumed  to  be  overdense 
or  to  be  operating  below  the  plasma  frequency.  The  thickness  of  the  sheath 
is  assumed  to  be  very  small  compared  with  a  wavelength,  thus  making  radio 
transmission  possible. 

In  1962  it  was  observed  (Ref.  1)  that  the  complex  dielectric  constant  of  a 

plasma  slab  could  be  simulated  by  a  rodded  medium.  Physically,  the  currents 

induced  in  the  rods  behave  like  the  oscillatory  currents  in  a  plasma  making 

the  simulation  possible.  In  many  cases,  however,  simulating  a  thin  over- 

dense  plasma  with  an  artificial  dielectric  is  difficult  because  of  fabrication 

limitations.  For  example,  the  simulation  of  a  sheath  with  a  thickness  equal 

to  1  /  5 0  requires  the  spacing  between  adjacent  rods  to  be  equal  to  or  less 

than  a  1  /  50  \  . 

o 

It  will  be  shown  that  a  wire  grating  structure  behaves  like  a  thm  overdense 
plasma  by  simulating  the  equivalent  surface  current  density  of  the  sheath.  It 
is  obvious  that  the  wires  cannot  simulate  all  the  complexities  of  an  actual 
re-entry  plasma,  but  at  this  stage  the  simulation  of  the  simplest  properties 
can  prove  to  be  very  useful  in  antenna  evaluation  and  design. 


II.  PLASMA  SURFACE  IMPEDANCE 


The  first  step  in  the  development  is  to  understand  that  which  characterizes 
the  reflection  and  transmission  properties  of  the  sheath.  For  the  purpose  of 
this  analysis,  the  sheath  is  assumed  to  be  an  infinite  slab  with  a  uniform  plane 
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wave  reflecting  off  its  surface  at  an  angle  0,  as  shown  in  Fig.  la.  When  the 
sheath  is  overdense  and  very  thin,  it  looks  like  an  inductive  surface  current 
to  the  incident  electromagnetic  signal.  Therefore,  it  should  be  possible  to 
simulate  the  sheath’s  characteristics  by  replacing  it  with  an  equivalent  sur¬ 
face  current  density.  It  is  interesting  to  note  that  J.  Wait  (Ref.  2)  and  other 
authors  have  used  this  technique  to  obtain  approximate  solutions  to  complex 
boundary  value  problems. 

The  jump  condition  used  by  Wait  can  be  developed  by  applying  Ampere’s  Law 
and  Faraday’s  Law  to  a  rectangular  circuit  as  shown  in  Fig.  lb.  The  appli¬ 
cation  of  Faraday’s  Law  leads  to  the  continuity  of  the  tangential  E  field  across 
the  sheath.  The  surface  current,  which  is  equal  to  the  discontinuity  of  the 
tangential  H  field,  is  obtained  by  using  Ampere's  Law. 

AEt  =  0  (1) 

AHt  =  Js  (2) 

The  surface  current  Jg  is  determined  by 

J  =  o  d  E  (3) 

si  '  ’ 

where  cr  is  the  complex  conductivity  of  the  sheath,  d  is  the  sheath  thickness, 
and  Ej  is  the  tangential  electric  field  at  the  surface.  The  reflection  and 
transmission  properties  of  the  sheath  are,  therefore,  governed  solely  by  the 
equivalent  surface  current  density  Jg.  The  surface  impedance  i-s  related 

to  Eq.  (3)  through  Ohm’s  Law 


Z 


(ohm  s/  square) 


-2- 


I 


Y 

Fig.  la.  Equivalent  Sheath  Problem  for  Jump  Conditions 
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Fig.  lb.  Slab  Problem 
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The  resulting  surface  impedance  for  a  thin  plasma  is 


Z 

sp 


trd 


(4) 


The  well-known  expression  for  conductivity  of  a  cold  plasma  (Ref.  3)  is 


a 


2 

€  CO 

°  P 

vc  +  jw 


where  is  the  collision  frequency  and  co^  is  the  plasma  frequency.  The 
surface  impedance  reduces  to 
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sp 


V  +  1  CO 

c  J 


2, 

€  CO  d 

o  p 


when  cr  is  substituted  into  Eq.  (4).  The  last  equation  can  be  expressed  in 
terms  of  three  nondimensional  quantities,  W,  V  and  D. 
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D  =  k  d 
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The  free  space  propagation  constant, 


k  ,  is  defined  as  follows 
o 


_  CO  _  2  7T 
o  c  X 

o 

The  surface  impedance  of  the  sheath  is  equivalent  to  a  series  RL  circuit. 

The  resistive  term  is  associated  with  the  loss  mechanism  in  the  plasma,  and 
the  reactive  part  is  related  to  the  phase  difference  between  the  motion  of  the 
electrons  and  the  applied  electric  field. 

The  basic  assumptions  associated  with  the  jump  conditions  (Ref.  4)  are  (1) 
the  frequency,  go,  and  the  thickness,  d,  must  be  such  that  the  flux  linking  the 
rectangular  circuit  is  negligible,  and  (2)  the  fields  inside  the  plasma  have  only 
z  dependence.  These  assumptions  are  satisfied  when 


GO 


and 


k  d  <<  1 
o 


HI.  WIRE  GRID  STRUCTURE 


The  behavior  of  a  plane  of  equally  spaced  wires  has  been  discussed  at  some 
length  by  various  authors  (Refs.  5,  6  and  7).  In  the  case  under  consideration, 
the  electric  field  is  assumed  to  be  parallel  to  the  wires  as  shown  in  Fig.  2. 
When  the  diameter  of  the  wires  is  much  less  than  the  wire  spacing  and  the 
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spacing  is  loss  than  a  half  a  wavelength,  then  the  surface  impedance  of  the 
grid  structure  (Refs.  8  and  9)  can  be  shown  to  equal 


R  +  jZ  *ln(-^-) 
s  J  o  K  \2rrr  / 


where  a  is  the  spacing  between  adjacent  wires,  r  is  the  radius  of  the  wires, 
and  is  the  free  space  wavelength.  If  the  wire  spacing  becomes  greater 
than  half  a  wavelength,  the  grid  structure  starts  to  take  on  the  character  of 
an  optical  grating  and  the  appearance  of  more  than  one  reflected  beam  is 
sometimes  observed.  The  internal  surface  resistance,  R  ,  of  a  grid  struc¬ 
ture,  consisting  of  n  wires  per  unit  width,  is  equivalent  to  the  resistance  per 
unit  length  of  n  wires  in  parallel.  It  follows  that 


R  -  aR 
s  o 


where  R^  is  the  internal  resistance  per  unit  length  of  one  wire. 


IV.  EQUIVALENCE  BETWEEN  PLASMA  AND  PERIODIC  STRUCTURE 


Next  in  the  development  is  the  comparison  of  the  surface  impedance  of  the 
wires  with  that  of  the  plasma  sheath.  The  surface  impedance  of  the  plasma 
sheath  is  given  by  Eq.  (5)  and  the  wire  grid  structure  by  Eq.  (9). 


Z 

sp 


Z  V  z 

-  °_  +  J _ 2_ 

2  ^  J  2 
W^D  W^D 


Z  =  R 
sw  s 


+  jZ  ^ 
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If  one  notes  that  both  impedances  have  a  resistance  and  an  inductive  reactance, 
then  the  relations  connecting  the  wire  grid  and  plasma  sheath  parameters  are 
found  by  equating  the  real  and  imaginary  parts  of  the  respective  impedances. 


R 
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Z 

o 


V 

w2d 


1 

w2d 


Equations  (11)  and  (12)  are  then  reduced  to  the  following: 


R 

s 
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C  V 

c 
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a  ln(z¥?)  2ttc2 

X  X  2, 

o  X  u,  d 

o  p 


(ID 


(12) 


1 

c  "  ~~ 

v/p  e 
v  o  o 

indicating  that  the  simulation  is  independent  of  frequency.  Equations  (11)  and 
(12)  constitute  the  equations  used  to  design  a  wire  grid  structure  having  an 
equivalent  W,  D,  and  V. 


V.  PRACTICAL  CONSIDERATIONS 


High  collision  frequencies  require  large  surface  resistance  as  seen  from 

Eq.  (11).  The  surface  resistance,  R  ,  must  be  of  the  order  of  500  ohms  per 

s 
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square,  if  collision  frequencies  in  the  neighborhood  of  the  operating  frequency 
are  desired.  One  promising  method  of  obtaining  large  internal  resistance 
consists  of  coating  quartz  fibers  with  a  thin  metallic  layer  of  thickness  t,  (see 
Fig.  3a).  By  making  the  thickness  of  the  metallic  layer  much  less  than  a  skin 
depth,  the  resulting  internal  resistance  is  independent  of  frequency  and  is 
given  by 


R  = - - - -  (ohms/meter)  (13) 

o  a  Zirrt 
m 

where  u  is  the  conductivity  of  the  metallic  layer.  With  this  technique, 
resistances  of  1000  ohms  per  inch  have  been  obtained. 

The  equally  spaced  wires  will  only  simulate  the  inductive  character  of  the 
plasma  when  the  electric  field  is  parallel  to  the  wires.  When  the  electric 
field  is  perpendicular  to  the  wires,  the  incident  wave  is  practically  unaltered 
by  their  presence.  A  more  detailed  analysis  (Ref.  10)  shows  that  the  surface 
admittance  is  small  for  this  polarization.  A  thin  plasma  sheath  can  be  simu¬ 
lated  by  a  square  lattice  screen  structure  independent  of  the  polarization 
(see  Fig.  3b).  The  incident  electric  field  can  be  separated  into  two  com¬ 
ponents,  each  one  parallel  to  one  set  of  wires.  The  resulting  surface  imped¬ 
ance  is  given  by  Fq.  (9). 

VI.  REFLECTION  AND  TRANSMISSION  COEFFICIENTS 


The  reflection  and  transmission  coefficients  for  the  plane  of  wires  and  the 
plasma  sheath  can  be  derived  from  the  jump  conditions  given  in  Eqs.  (1)  and 
(Z).  Assume  that  a  plane  wave  is  incident  at  an  angle  0  from  the  left  and  the 
electric  field  is  parallel  to  the  plane  of  the  interface  as  shown  in  Fig.  Z. 
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Application  of  the  jump  conditions,  Eqs.  (1)  and  (2),  to  this  geometry  yields 


(H1  +  Hr 


0  T 


El  +  Er 


(14) 


( E1  +  Er  -  El)  =  0 


(15) 


Since  the  incident,  reflected,  and  transmitted  fields  are  uniform  plane  waves, 
the  electric  field  is  related  to  magnetic  field  through  an  intrinsic  impedance 
Use  of  this  concept  reduces  Eq.  (14)  to 


(E1  -  Er 


_t.  cos  O  E1  +  Er 

E  >  ~2—  =  — - 

o  s 


(16) 


Solution  of  Eqs.  (15)  and  (16)  for  the  ratio  Er/EX  leads  to  a  reflection 
coefficient  given  by 


r 


1  + 


-  1 

2  Z  cost) 
s 

Z 

o 


(17) 


In  a  similar  manner,  the  transmission  coefficient  is  obtained  by  solving  Eqs. 
(15)  and  (16)  for  the  ratio  eVe1 


t  = 


1  4 


2Z  cos  0 
s 


(18) 


The  reflection  and  transmission  coefficients  for  the  thin  plasma  sheath  are 
obtained  by  substituting  the  surface  impedance  given  by  Eq.  (5)  into  Eqs.  (17) 
and  (  18). 
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(19) 


P  i  ±  2V  cos  0  2  cos  O 

1+  - 5 -  +J  - 2 

W  D  W  D 

2V  cos  0  2  cos  O 

2  +  J  2 

,  =  W^D _ \v  D 

p  t  i  2V  cos  O  .  2  cos  S 

F  1+  - 1 -  +  J  - 2 - 

W  D  W  D 


(20) 


Similarly,  the  reflection  and  transmission  coefficients  for  the  wire  grid 
structure  are  given  by 


!  +  ^  cos  0  +  j  ^  cos  Oln^) 
o  o 

7^  cos  O  +  j^cos  Gln^) 

o _ o_ _ 

1  +  cos  0  +  j  IT  cos  0  ln(l5?) 

o  o 


(21) 


(22) 


The  subscripts  p  and  w  denote  the  quantities  pertaining  to  the  thin  plasma 
sheath  and  the  wire  grid  structure,  respectively. 


VII.  EXPERIMENTAL  DATA 


The  reflection  and  transmission  properties  of  different  wire  grid  structures 
were  measured  in  an  X-band  waveguide  geometry.  The  dominant  TE^  mode 
in  a  rectangular  waveguide  is  analogous  to  two  uniform  plane  waves  bouncing 
off  the  walls  of  the  waveguide  as  shown  in  Fig.  4  (Ref.  11).  The  equivalent 
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Fig.  4*  Grid  Structure  in  X-band  Waveguide 
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angle  of  incidence,  O,  of  these  two  waves  is  related  to  the  waveguide 
impedance  and  waveguide  wavelength  as  shown  below: 


Z 

g 


Z 

o 

cos  O 


(23) 


X 

g 


X 

o 

cos  o 


(24) 


With  this  in  mind,  the  measurements  obtained  in  the  waveguide  geometry  are 
equivalent  to  the  same  measurements  performed  in  free  space  but  at  an  angle 
of  incidence  given  by  Eqs.  (23)  and  (24).  The  equations  governing  the  complex 
reflection  and  transmission  coefficients  for  the  wire  structure  in  the  wave¬ 
guide  geometry  are  obtained  by  substituting  Eqs.  (23)  and  (24)  into  Eqs.  (21) 
and  (22). 
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w 


1  +  2(Rs/Zg)  +  jha  ln(^:)/(Xg) 


2(Rs/Zg)  +  j 

2a  ln| 

^r) 

|/(Xg) 

1  +  2(Rs/Zg)  +  j 

2a 

H 

life) 

|/(Xg) 

(25) 


(26) 


The  experimental  data  were  measured  with  an  X-band  waveguide  interfer¬ 
ometer  shown  in  Fig.  5.  The  readings  on  the  precision  attenuator  and  phase 
shifter  are  related  to  the  complex  transmission  coefficient  of  the  sample 
being  tested. 

The  slotted  line  located  directly  in  front  of  the  test  section  yields  information 
pertaining  to  the  reflection  coefficient.  The  position  of  the  nulls  of  the  stand¬ 
ing  wave  pattern  are  related  to  the  reflection  phase,  and  the  VSWR  reading 
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CRYSTAL 
detector 
(  H PX  485B 


yields  the  magnitude  of  the  reflection  coefficient.  When  the  phase  shift 
measurements  are  being  performed,  the  probe  antenna  should  be  extracted 
from  the  slotted  line. 

The  actual  experiment  consisted  of  measuring  both  the  complex  reflection 
and  transmission  coefficients  for  wire  grid  structures  having  a,  ranging 
from  0.  15  in.  to  0.  3  in.  and  r,  ranging  from  0.  75  to  4.  0  mils.  The  wave¬ 
guide  test  section  which  holds  the  wire  grid  structure  in  place  is  shown  in 
Fig.  6.  The  lossless  data  (Rg  =  0)  are  plotted  as  a  function  of 

\g  W2P 

-  Ms?) ' cos  ° 


and  compared  with  Eq.  (25).  The  equivalent  plasma  frequency  and  thickness 
of  the  simulated  sheath  can  be  obtained  from  Eq.  (  1  2).  Figures  7  and  8 
summarize  the  measurements  performed  on  the  lossless  copper  wire  models. 
As  the  diameter  of  the  wire  increases  or  the  spacing  decreases,  the  reflec¬ 
tions  off  the  grid  structure  become  more  severe  as  is  seen  in  Fig.  7.  The 
simulation  of  more  dense  plasma  layers  can  therefore  be  accomplished  either 
by  using  larger  wires  or  by  increasing  the  number  of  wires  in  the  grid 
structure. 

The  experimental  data  for  the  lossy  coated  fibers  are  compared  with  Eqs.  (25) 
and  (26)  in  Figs.  9  through  16.  The  lossy  wires  are  quartz  fibers  ranging 
from  1.  5  to  4.  0  mils  in  diameter  with  a  gold  pottery  glaze  fired  on  their 
surface.  The  value  of  R  is  obtained  by  measuring  the  linear  resistance,  R  , 
with  a  dc  ohmmeter  and  relating  Rq  to  Rg  by  means  of  Eq.(io).  The  addition 
of  losses  in  the  grid  structure  tends  to  decrease  the  severity  of  the  reflec¬ 
tions,  thereby  increasing  the  transmission  of  electromagnetic  energy  as 
shown  by  the  experimental  data.  The  decrease  in  the  reflection  coefficient  is 
attributed  to  the  improved  matching  between  the  simulated  sheath  and  its 
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Fig.  7.  Lossless  Reflection  Magnitude  Data 
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Transmission  Phase  Data  for  an  Equivalent  Collision 
F  requency  of  1.476  Gc 


Collision  Frequency  of  12.82  Gc 
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Transmission  Phase  Data  for  an  Equivalent  Collision 
Frequency  of  126  Gc 


surroundings.  As  the  losses  are  increased,  the  surface  impedance  of  the 
simulated  sheath  increases  which,  in  turn,  makes  the  reflections  smaller. 

Additional  experimental  data  were  obtained  for  a  wire  mesh  structure.  How¬ 
ever,  in  this  case,  the  wires  making  up  the  mesh  were  not  orthogonal  with 
the  waveguide  walls.  The  results  indicate  no  measurable  deviation  in  the 
reflection  characteristics  from  the  models  having  wires  parallel  to  the 
electric  field.  The  reflection  and  transmission  coefficients,  therefore,  obey 
the  relationships  given  by  Eqs.  (25)  and  (26)  and  are  independent  of  the 
polarization  of  the  incident  fields. 


VIII.  SUMMARY 


The  simulation  of  a  thin  plasma  sheath  can  be  accomplished 
lattice  screen  structure  when  the  following  assumptions  are 


by  using  a  square 
valid 


CJ 


U> 


k  d  <  c  1 
o 

Large  collision  frequencies  can  be  simulated  by  using  coated  quartz  fibers. 

Collision  frequencies  in  the  neighborhood  of  the  operating  frequency  are 

realizable  by  this  method.  The  experimental  data  confirm  that  the  grid 

structure  will  simulate  the  impedance  characteristics  of  a  thin  plasma  sheath. 

2 

The  values  of  W  D  and  V  for  the  coated  fiber  grid  structures  are  in  a  range 
which  characterizes  the  sheaths  developed  around  slender  body  re-entry 
vehicles.  Plasma  sheath  simulation  with  a  grid  structure  could  therefore  be 
useful  in  ground  testing  and  evaluation  of  re-entry  antenna  systems. 
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